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ABSTRACT 



A deep mid- and far-infrared survey in the Extended Groth Strip (EGS) area 
gives 3.6 to 8/im flux densities or upper limits for 253 Lyman Break Galaxies 
(LBGs). The LBGs are a diverse population but with properties correlated with 
luminosity. The LBGs show a factor of 30 range in indicated stellar mass and 
a factor of 10 in apparent dust content relative to stellar mass. About 5% of 
LBGs are luminous at all wavelengths with powerful emission at rest 6/im. In 
the rest 0.9 to 2/zm spectral range these galaxies have stellar spectral slopes with 
no sign of an AGN power law component, suggesting that their emission is mainly 
powered by intensive star formation. Galaxies in this luminous population share 
the infrared properties of cold SCUBA sources: both are massive and dusty 
starburst galaxies at 2 < z < 3; their stellar mass is larger than lO n M . We 
suggest that these galaxies are the progenitors of present-day giant elliptical 
galaxies, with a substantial fraction of their stars already formed at z ~ 3. 
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1. Introduction 

The two most efficient ways of identifying galaxies at high redshifts are first, the Lyman 
break broadband-droput technique, and second, observations at submillimeter wavelengths. 
For example, surveys with the Submillimeter Common-User Bolometer Array (SCUBA) on 
the James Clerk Maxwell Telescope reveal high-redshift objects via emission from the cold 
dust they contain (Ivison et al. 2002; Smail et al. 2002; Chapman et al. 2003). Conventional 
wisdom is that the submillimeter observations preferentially select starburst galaxies at z > 2 
because the far-infrared luminosity peak shifts into the submillimeter band (Ivison et al. 
2002; Smail et al. 2002; Chapman et al. 2003). Most SCUBA sources are too faint in the 
optical and near-infrared bands for spectroscopic identification, but Chapman et al. (2003) 
and Simpson et al. (2003) have confirmed that most SCUBA sources in their sample are at 
z > 2. Despite all the effort that has gone into SCUBA searches, Lyman Break Galaxies 
(LBGs) still constitute by far the largest well-studied rest-frame-UV selected galaxy samples 
at high redshift (Steidel et al. 2003). 

The relation between LBGs and SCUBA galaxies remains unclear, but recently Chap- 
man et al. (2005) have confirmed that some SCUBA galaxies have the typical rest-frame 
UV colors as LBGs at z ~ 2 (Steidel et al. 2004). Understanding both populations requires 
a direct comparison of properties, such as stellar mass and dust emission. The four bands 
viewed by the Infrared Array Camera (IRAC) aboard Spitzer, from 3.6^m to 8.0/xm, cover 
spectral energy distributions (SEDs) for galaxies at 2 < z < 3 in the rest-frame near-infrared, 
where the luminosities are good indicators of stellar mass (Cole et al. 2001; Drory et al. 2004; 
Glazebrook et al. 2004; Bundy et al. 2005). The Multiband Imaging Photometer for SIRTF 
(MIPS) 24/zm band is an ideal probe of the PAH and hot dust emission from LBGs and 
SCUBA galaxies. Preliminary studies show that Spitzer can easily detect both LBGs and 
SCUBA galaxies (Barmby et al. 2004; Egami et al. 2004; Huang et al. 2004). 

This Letter presents the first study of a large LBG sample in the Extended Groth 
Strip (EGS) region based on Spitzer observations carried out with deep IRAC and MIPS 
photometry. We report discovering a new type of LBG, the Infrared-Luminous LBG (ILLBG) 
with strong infrared emission at rest-frame 6//m, whose properties are very similar to those 
of SCUBA galaxies. ILLBGs are massive, dusty, and powered mainly by extremely rapid 
star formation. 
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2. Spitzer Observations 

IRAC (Fazio et al. 2004a,b) observations of the EGS were carried out in two epochs: 
2004 January and June, covering a 2° x 10' strip. The MIPS observations were done in 2004 
June in scanning mode. Limiting flux densities (5<r) are roughly 0.5yuJy at 3.6 and 4.5/im, 
2.7//Jy at 5.8 and 8.0^m, and 60/iJy at 24/zm. We use the same data reduction procedure 
for both IRAC and MIPS data as Huang et al. (2004) and Egami et al. (2004). 

The LBG sample in the EGS area was taken from Steidel et al. (2003). Among a total 
of 334 objects in the Steidel et al. (2003) catalog, 188 are spectroscopically identified as 
galaxies at z~3, 3 are QSOs, 3 are AGN, and 7 are stars. The remaining 133 objects are not 
identified. In the photometric catalog 253 objects are in the area covered by both IRAC and 
MIPS imaging. Of these, 211 are detected at 3.6/xm, 199 at 4.5/um, 53 at 5.8/xm, and 44 at 
8.0/xm. We detect 11 LBG counterparts in the MIPS 24^m image including all 3 QSOs and 1 
AGN identified by Steidel et al. (2003). 1 The 24/xm-detected counterparts are also detected 
at all shorter wavelengths and are among the brightest sources in the sample. We define the 
LBGs with 24 fim detections (/24^m > QOfiJy) as Infrared-Luminous Lyman Break Galaxies 
(ILLBGs). The 24/im limiting flux is equivalent to L^^m = 6 x 10 8 L Q ^imT 1 for galaxies at 
z = 3, or L ir = 3 x 10 10 L© if we use the M82 SED to convert the luminosity at 6.2/xm to 
the total infrared luminosity. 

We compare the LBG sample with a set of Lockman Hole area SCUBA sources having 
secure identifications (Huang et al. 2004; Egami et al. 2004). 2 Egami et al. (2004) classified 
these SCUBA sources into two types (cold and warm) based on the dust temperature inferred 
from their SEDs. The cold sources exhibit the 1.6yum H~ opacity minimum bump in the 
IRAC bands, a stellar feature typically seen in galaxy SEDs. The cold SCUBA sources are 
well approximated by the SED of Arp 220, a cold dusty starburst. The warm sources are 
AGN-dominated with a power-law continuum SED, similar to that of Mrk 231. There are 
four cold and two warm sources in the Lockman Hole SCUBA sample (Egami et al. 2004; 
Huang et al. 2004). 



1 Two additional LBGs may be detected at 24/xm, but there are non-LBG galaxies within 3" of them on 
the sky. Given the angular resolution of MIPS, we cannot be sure which objects are responsible for the 24^m 
flux, and to be conservative we do not count these two as 24^m LBG detections. 

2 Finding counterparts to SCUBA sources is difficult because there are many faint optical and infrared 
sources within the SCUBA error circle. For reliable identifications, radio observations are needed to give 
accurate positions. Four SCUBA sources in the EGS IRAC+MIPS area have good radio observations (?), 
but none of the counterparts thereby identified has a known rcdshift. In fact, SEDs suggest that at least two 
of the four are at z < 2. Thus the EGS SCUBA sources do not provide a useful comparison sample, and we 
compare the EGS LBGs with the Lockman Hole SCUBA sources instead. 
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3. Infrared Colors of LBGs 

The LBGs exhibit a much wider range of flux densities in the IRAC bands than in the K 
band. Figure 1 shows that the range of flux densities for the LBGs spans 4 magnitudes 

compared to only 1.5 magnitudes for the range of K-band flux density seen by Shapley et 
al. (2001). Figure 1 also shows a correlation between [3.6] and R — [3.6] for LBGs. This 
relation is partially due to a selection effect (all LBGs are brighter than R = 25.5 (Steidel 
et al. 2003)), but there is an obvious absence of LBGs with blue colors at higher 3.6yum flux 
densities. All ILLBGs (as defined above) are redder than R — [3.6] = 1.5, and most have 
R — [3.6] > 2.4. Three ILLBGs are so red that they are close to qualifying as extremely red 
objects (ERO) (R - [3.6] > 4, Wilson et al. 2004). LBGs detected at 8/iin are redder than 
other galaxies {(R — [3.6]) = 1.95 versus 0.94 for LBGs detected at 3.6/um but not 8yum) but 
not so red as the ILLBGs. All six Lockman Hole SCUBA sources have red R — [3.6] colors 
(Wilson et al. 2004), comparable to those of the ILLBGs. 

Figure 2 compares the observed LBG colors with those predicted by two simple stellar 
population synthesis models (single burst and constant star formation). These two models 
correspond to two extreme cases of a model having a exponentially decay star formation 
rate, ip(t) = e~*/ r , where r = is the single burst model and r = oo is the constant star 
formation model. For any given E(B — V), the constant star formation model predicts 
the reddest R — [3.6] color for galaxies at z ~ 3. A combination of the two models with 
varying amounts of dust attenuation can reproduce almost all the LBG colors. Shapley et 
al. (2001), Papovich et al. (2001), and Labbe et al. (2005) argue that a model with constant 
star formation rate fits the SEDs of near-infrared bright LBGs better. This model should 
also be very close to the real star formation history for the 8/j.m bright LBGs and ILLBGs. 
As shown in Figure 2, g' — R is more sensitive to dust extinction, and R — [3.6] is sensitive 
to both age and dust extinction for galaxies with constant star formation history. A rough 
estimate of the extinction range as 0.1 < E(B — V) < 0.4 is consistent with what Shapley et 
al. (2001) and Papovich et al. (2001) measured in their LBG samples. The LBGs with 8/im 
detections tend to be older but little dustier than those not detected at 8/xm. The ILLBGs 
are more extreme, most very red in g' — R and R — [3.6], some beyond the model prediction. 
Shapley et al. (2001) also found that the SEDs of LBGs with the reddest R — K s defied 
any simple model fitting and argued that those are starburst galaxies with multi-component 
morphologies, each with different colors. Such a red color could also be due to a significant 
contribution from an obscured AGN, but we argue below that significant AGN emission is 
inconsistent with the [4.5] — [8.0] colors. 
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4. Origin of 24/xm emission from ILLBGs 

The observed 24^m emission («6/xm rest wavelength) from the ILLBGs can come from 
either the power law component of an AGN (e.g., Elvis et al. 1994; Egami et al. 2004; Stern et 
al. 2005) or from warm dust heated by a starburst (e.g., Sanders & Mirabel 1996; Rigopoulou 
et al. 1999). For strong starbursts, the 6.2/xm PAH feature can enhance the observed flux 
by perhaps a factor of 1.5 (Tran et al. 2001; Peeters et al. 2004). As Ivison et al. (2004) 
point out, a key diagnostic for distinguishing between starbursts and AGN is the existence 
of a change of slope around rest 3-4/iin: strong AGN have nearly a constant slope from 2 to 
10/im, whereas starbursts show a distinct minimum in the 3-4/im range. Translated to an 
observational criterion, rest J — K s > 2 implies an AGN (Cutri et al. 2001; cf. Ivison et al. 
2004 Fig. 3). 

Figure 3 shows the [8.0] — [24] vs J—K s color-color diagram for starburst- AGN separation 
suggested by Ivison et al. (2004). QSOs and warm SCUBA sources exhibit rest-frame J — K s 
color as red as local AGNs (Cutri et al. 2001), while all ILLBGs except one have J —K s < 1.5. 
Therefore we suggest that the 24/xm flux densities of ILLBGs are mainly due to strong dust 
emission caused by intensive star formation. ILLBGs have much redder [8.0] — [24] colors than 
QSOs and AGNs, steeper than typical power-law SEDs but consistent with strong starbursts. 
The SCUBA sources occupy the same overall portion of color-color space in this plot as do 
ILLBGs, AGNs, and QSOs at z ~ 3. In other words, the cold SCUBA sources appear to be 
Arp-220-like starburst galaxies at the upper-left corners, and the warm SCUBA sources are 
AGNs. We conclude that ILLBGs, like cold SCUBA sources, are starburst galaxies. 

5. Characterizing the Galaxy Population at z = 3 

The faint blue LBGs that lack 24yum detections could have dust emission with flux 
densities below our detection threshold. For example, even blue dwarf galaxies are found 
to have dust emission at rest-frame 6//m, presumably with PAH features (Rosenberg et al. 
2005). We used a stacking technique to assess the average dust emission for the blue LBG 
population at z=3. The LBGs without 24yum detection were divided into two groups for 
stacking: those with 8^m detections and those without. There are a total of 198 LBGs 
lacking 24/xm detections including both spectroscopically identified and unidentified ones. 
In this subset, 31 LBGs have 8yum detections. The 8/xm and 24/im images for both groups 
were stacked separately. 

The stacking results show a significant difference between LBGs with and without 8yum 
detections. The LBGs without 8/xm detections have (Mft-(Uega)) = —21.5, right in the range 
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exhibited by local dwarf galaxies as shown in Figure 4. The upper limit in the [8.0] — [24] 
colors is also consistent with colors of local blue dwarf galaxies. This strongly suggests that 
most faint blue LBGs in the sample have similar stellar masses, dust obscuration, and star- 
formation histories as dwarf galaxies in the local universe. The LBGs with 8/xm detections, 
including the ILLBGs, have {MxiVega)) = —24.5, brighter than the local M^iVega) = 
-23.4+5 *log(/i) (Cole et al. 2001; Huang et al. 2003). The (M K (Vega)) and ([8.0] -[24]) for 
the 8/xm-detected LBGs are only slightly fainter and bluer than those of ILLBGs, suggesting 
that they may be similar kinds of galaxies but a little fainter and bluer, just below the 24/xm 
limiting flux density. 

All the results show that the LBG population is far from uniform. Figures 1 and 4 
demonstrate that the LBGs with the largest stellar luminosities, as indicated by 8/xm flux 
densities, are older and dustier than the LBGs with lower stellar luminosities. A larger LBG 
sample with infrared photometry and deeper 8^m observations are needed to tell whether 
there are two distinct populations or a continuous range of galaxy properties. 

Rigopoulou et al. (2005) conducted a detailed analysis of the stellar populations and 
estimated stellar masses for the 8/xm-selected LBGs. All were found to be massive galaxies(> 
4 x 10 10 M o ) with the ILLBGs tending to have higher masses (> lO n M ). These masses are 
consistent with the mean dynamic mass of 1.5 x 10 n M Q for starburst-type SCUBA galaxies 
(Swinbank et al. 2004; Tecza et al. 2004). Massive galaxies with 1O 11 M have also been found 
at 1 < z < 2 (Cimatti et al. 2002; Daddi et al. 2004; Saracco et al. 2004; Glazebrook et al. 
2004; Tecza et al. 2004; Labbe et al. 2005). These observations challenge current theoretical 
predictions of galaxy formation (Kauffmann et al. 1999; Shu et al. 2001; Baugh et al. 2003; 
Somerville et al. 2004). This study provides further confirmation that at least some massive 
galaxies were formed by z = 3, favoring galaxy formation and evolution models assuming 
more rapid stellar mass accumulation scenarios (Nagamine et al 2005). 



6. CONCLUSIONS 

This infrared study of LBGs reveals that they are a diverse group in terms of both their 
masses and dust contents. The observed 8/iin (rest 2/im) flux densities, which are roughly 
proportional to stellar mass, extend over at least 1.5 orders of magnitude, and the ratio of 
24 to 8^m flux densities, which measure the fraction of ISM emission, extends over more 
than an order of magnitude. 

Among 253 LBGs in the EGS area covered by Spitzer, 5% were detected with /W^m > 
60/iJy. We refer to these as ILLBG. ILLBGs have much redder observed R— [3.6] colors than 
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the other LBGs in the sample. The rest-frame J — K s colors for most ILLBGs are stellar 
with no sign of an AGN power-law component. Therefore we suggest that the 24^m emission 
is dominated by dust heated in a massive starburst, probably including some contribution 
from the PAH feature at rest-frame 6.2/im. ILLBGs and cold SCUBA sources share the 
same infrared properties, suggesting that they are closely related and may belong to the same 
population. If so, it ought to be possible to detect submillimeter emission from ILLBGs. Both 
ILLBGs and cold SCUBA sources tend to be massive with a typical stellar mass of 1O 11 M . 
Such masses make them candidates for the progenitors of present-day giant ellipticals. 

In contrast to ILLBGs, most LBGs in the sample are only detected at 3.6 and 4.5yum 
and have faint mid-infrared luminosities and blue R — [3.6] colors. Their lower average 8/xm 
flux density measured with the stacked image implies that the faint LBGs have much less 
stellar mass than ILLBGs, in a range similar to local dwarf galaxies. An upper limit of 
[8.0] — [24] color for average faint LBGs is also consistent with dust emission of local dwarf 
galaxies. 
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lowship under grant AST-0302049. This work is based on observations made with the Spitzer 
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20 21 22 23 24 25 26 

[3.6] AB 

Fig. 1. — Observed-frame R — [3.6] versus [3.6] color-magnitude diagram for all LBGs in 
our EGS sample. Black circles indicate LBGs detected only at 3.6/iin and 4.5/im. The 
green circles are LBGs detected at 8/xm but not 24/zm, and red circles are LBGs detected 
in both 8//m and 24/xm bands. For all these objects, filled symbols designate sources having 
spectroscopic redshifts while open symbols are for objects with photometric redshifts only. 
Filled squares and diamonds represent the cold and warm SCUBA sources in the Lockman 
Hole (Egami et al. 2004). The straight line is the R-band limiting magnitudes of R=25.5 
(Steidel et al. 2003). 
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Fig. 2. — Observed-frame g' — R vs R — [3.6] color-color diagram for all LBGs detected 
with IRAC. Symbols are as in Figure 1. The lines show colors from stellar synthesis models: 
dashed lines are the single burst models, and solid lines are the constant star formation 
models. The blue lines are dust-free models, and the red lines are models with E(B—V) = 0.3 
and a Calzetti et al. (2000) reddening law. The models run from 0.1 Myr at the blue end of 
the lines for the reddened models, off the plot for the dust-free models to 2 Gyr as indicated 
for the constant star formation models, above the top of the plot range for the single burst 
models. Maximum ages for the single burst models are indicated near the top of the figure. 
Overall the set of models spans nearly the entire color range of the data. 
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Fig. 3. — Observed-frame [8.0] — [24] vs rest-frame J — K s color-color diagram for ILLBGs, 
AGNs, and SCUBA sources. At z=3, observed [4.5] — [8.0] is nearly equivalent to rest-frame 
J — K s , and we have converted to J — K s using K-corrections from Hogg et al. (2002). The 
J — K s color is in Vega magnitudes for comparison with the 2MASS results (Cutri et al. 
2001; Cole et al. 2001), but [8.0] - [24] is in AB magnitudes. 
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Fig. 4. — Observed-frame [8.0] — [24] vs rest-frame MxiVega) color-magnitude diagram 
for ILLBGs, AGNs, and SCUBA sources. We derived the observed-8/imto rest-K-band K- 
correction with a spiral SED for ILLBGs and cold SCUBA sources and an AGN SED for 
AGNs, QSOs, and warm SCUBA sources. A mean redshift of z = 2.97 for LBGs was assigned 
to those LBGs without spectroscopic redshifts. The upper limit shown at Mic(Vega)— -21.5 
is for the stacked LBGs lacking 8/iin detections. The green filled circle is for the stacked 
LBGs having 8/im detections. [8.0] — [24] color for objects at z ~ 3 is equivalent to the rest- 
frame color K — [6.0]. We plot the K — [5.8] color for local blue dwarf galaxies (Rosenberg et 
al. 2005) for comparison. The solid line is the MIPS 24/im limiting flux density for galaxies 
at z = 3. 



